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Abstract: We present calculations of the single-particle exci- 
^3 ■ tation spectrum for a 2D strong-coupling superconductor in a 

conserving approximation. Spectral weight at low frequency is 
q ! substantially reduced as the superconducting transition is ap- 

proached from the normal state. The suppression of low-energy 
excitations is a consequence of Cooper-pair fluctuations that 
are described self-consistently in the fluctuation exchange ap- 
proximation. The static and uniform electromagnetic response 
provides a measure of the super fluid density and a fully self- 
consistent indication of the superconducting transition temper- 
ature. 

The pseudogap observed in the normal state of the high-temperature superconductors 
has focused interest on fluctuations directly associated with superconductivity (Emery and 
Kivelson 1995, Randeria, et al. 1992, Doniach and Inui 1990). The occurrence of supercon- 
ductivity in two-dimensional copper-oxide layers encourages speculation that Cooper-pair 
fluctuations affect the low-energy excitation spectrum irrespective of the pairing mechanism 
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or of the symmetry of the order parameter. It has been proposed that fluctuations of the 
phase of the order parameter suppress the physical transition temperature, i.e. the tempera- 
ture at which macroscopic signatures of phase coherence first appear, below the temperature 
at which a gap forms in the single-particle excitation spectrum (Emery and Kivelson 1995, 
Trivedi and Randeria 1995). Here we investigate the excitation spectrum in a self-consistent 
calculation that includes fluctuations of both the amplitude and phase of the order parame- 
ter. 

The short coherence lengths characteristic of the high-temperature superconductors sug- 
gest that an appropriate description lies somewhere between the BCS picture and Bose 
condensation of bound fermion pairs (Randeria et al. 1989). The experimental observation 
of what is apparently a well defined Fermi surface in the high-temperature superconductors 
and the presence of a sharp drop in n(k) for \U\ ~ W/2 in quantum Monte Carlo calculations 
on the Hubbard model with an attractive two-body interaction (Trivedi and Randeria 1995), 
indicate that a pairing interaction of intermediate coupling strength on the scale of a typical 
electronic bandwidth may be of physical interest .[] Recent angle-resolved photoemission ex- 
periments on BSCCO have been interpreted as being consistent with a large Fermi surface 
and therefore a high electronic density (Ding et al. 1997). The nature of the crossover from 
the overlapping Cooper pairs of BCS theory to a Bose condensation of preformed pairs has 
been explored by several authors in the context of continuum and lattice models (Leggett 
1980, Nozieres and Schmitt-Rink 1985, Eagles 1969, Randeria et al. 1992, Luo and Bick- 
ers 1993, Micnas et al. 1990, and Sa de Melo et al. 1993). Taken together, these works 
highlight the importance of collective behavior in the particle-particle channel, which is ex- 
pected to play a significant role, particularly for large interaction strength, in determining 
the transition temperature, the nature of the phase transition, and, of central interest here, 
the electronic excitation spectrum. 

In a previous work (Deisz et al. 1998a), hereafter called DHS, we explored superconduc- 
tivity in the fluctuation exchange approximation (FEA) (Bickers et al. 1989), a conserving 
approximation (Baym 1962) beyond mean field theory. The superfluid density, its tempera- 
ture derivative, and the specific heat all show dramatic effects of fluctuations, and scale slowly 
with increasing lattice size. This slow scaling so far precludes an unambiguous determina- 
tion of the type of phase transition, and could even be consistent with a Kosterlitz-Thouless 
transition, but in any case it is of interest to know how Cooper-pair fluctuations affect 
the electronic excitation spectrum in a theory in which fluctuations are powerful enough to 
correctly renormalize the transition temperature from mean-field theory. To address this 
question, we calculate single-particle spectral functions in the FEA and show the tempera- 
ture evolution of the excitation spectrum side-by- side with the superfluid density calculated 
on lattices of 128 x 128 momentum points. We consider a density n = 0.75 roughly consistent 
with that obtained in electronic structure calculations. This provides a connection with the 
discussion of the nature of the transition presented in DHS which also shows that transition 
temperatures obtained on small lattices at this density are in good agreement with quantum 
Monte Carlo calculations. While it is expected that phase fluctuations are most important 
for the suppression of long-range order, the FEA contains fluctuations of both the amplitude 

1 By contrast, for a sufficiently strong pairing interaction, all electrons are expected to participate in the 
formation of 'diclcctronic molecules' and no Fermi surface is expected. 
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and phase of the order parameter. 

The central results of this paper are contained in Fig. 1 which shows the electronic spec- 
tral weight accessible to a photoemission experiment at four temperatures. Spectral weight is 
indicated by color in an energy (y-axis) momentum (x-axis) plane. The momenta shown are 
along the (1,0) direction from r to X in the square Brillouin zone. Red indicates the largest 
spectral weight; spectral weight decreases from red descending in the order of the colors of 
the spectrum to the smallest indicated by indigo. Only for the lowest temperature (Panel 1) 
is there a non-zero superfluid density; it is very small compared to that for zero temperature. 
The spectrum for the highest temperature, shown in Panel 4, is consistent with a band of 
short-lifetime excitations broadened by correlations and by thermal fluctuations. As the 
temperature is lowered, narrow yellow and red regions appear in the electronic excitation 
spectrum for some momenta above and below the Fermi surface. The 'sharpening' of these 
single-particle excitations^] signals the evolution and emergence of quasiparticle excitations^ 
with energies above and below the chemical potential \x. For momenta nearer to the Fermi 
surface, excitations do not 'sharpen' and the light green color indicates that spectral weight 
at low-energy decreases as the temperature is lowered. As we explicitly show below for a 
point very near the Fermi surface in the (1, 0) direction, spectral weight at zero energy is 
suppressed significantly as the transition is approached; relative to the 'sharp excitations' 
the maxima of the spectral functions associated with low-energy excitations are reduced by 
approximately 50%. Our results suggest a simple model discussed below. 

We begin with a model of coupled Cooper-pair fluctuations and electrons that consists 
of the Hubbard Hamiltonian with an attractive interaction U < threaded by a flux 

H((f> B ) = E l eX P("Tr) 4,a C r+x,<x + h.c.) 

r,a Vo 1 ^ 

~t Yl [ C l,aCr+y,a + h.c] ™r,T™r,|, (1) 

r,cr r 

together with the FEA. Here —t is the nearest-neighbor hopping matrix element, U(< 0) is 
an on-site attractive interaction and <p = hc/e. 

To self-consistently identify a superconducting phase transition in this model, we calcu- 
lated the internal energy E, and free energy F, as a function of an applied flux using the 
self-consistently determined self-energy and fully renormalized propagator. A finite super- 
fluid density D s leads to a non-zero curvature in E(<Pb) and F((/)b) at <pB = (Yang 1962, 
Fisher et al. 1973, Scalapino et al. 1992), 

Ffo) = F(0) + ^ s (T)A 2 + .., (2) 

E(<p B ) = m + \[D s -T^](^-r+... . (3) 

2 To be more precise, by 'sharpening' we mean that the peak in the spectral functions for these momenta 
becomes narrower and spectral weight at the maximum increases. 

3 While these excitations are relatively 'sharp,' the temperature dependence of their lifetime is not con- 
sistent with expectations for a Fermi liquid. This may not be surprising given the high temperature and the 
relatively large energy of these excitations. In any event, this is not the central point of this paper, and will 
not be pursued further here. 
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Figure 1: Single-particle spectral weight represented by color with darkest blue (indigo) corresponding 
to zero and red to ~ 1.4. Spectral weight increases in the order of the colors of the rainbow. Panel 1 is for 
T = 0.12; temperature increases by 0.02 going from panel to panel from left to right. The cc-axis in each 
panel indicates momentum from T- to X-points along the (1,0) direction and the y-axis indicates energy 
from —5.5 to 2.0. The white line indicates the position of the chemical potential. Spectral weight at low 
energy above the transition is suppressed with decreasing temperature (from 4 to 1) as indicated by the light 
green region sandwiched between the evolving red peaks. Each panel corresponds to an arrow in Fig. 2. 
Here, as in all figures, U = — 4 and the density is 0.75. 

This procedure is equivalent to obtaining D S (T) from the long- wavelength and zero-frequency 
limit of the fully self- consistent electromagnetic response function. DHS examined signatures 
of the phase transition in D S (T) and dD s (T)/dT for lattice sizes from 4 x 4 up to 64 x 64 
and for temperatures as low as 87 K (taking the bandwidth to be 1 eV). 

To study the effect of Cooper-pair fluctuations on the single-particle excitation spectrum, 
we examine the fully renormalized propagator and self-energy obtained from numerical solu- 
tions of the equations for the FEA with <pB = 0. A detailed description of the FEA together 
with a summary of our computational methods can be found elsewhere (Serene and Hess 
1991, Deisz, et al. 1994, Deisz, et al. 1998a, Deisz, et al. 1998b). Central to our discussion of 
the single-particle excitation spectrum is Dyson's equation, which relates the fully renormal- 
ized propagator G to the self-energy E and to the dispersion relation of the noninteracting 
system, 



where £k = ^k - 1 1 , A* is the chemical potential and e n = (2n+l)nT are Matsubara frequencies. 




(4) 



4 



The FEA self-energy includes, in addition to the Hartree term, the second-order term and 
the exchange of spin-fluctuations, density-fluctuations and Cooper-pair fluctuations, 

£(r, r) = U 2 [ Xph (r, r) + T pp (r, r) + T sf (r, r)} G(r, r) + U% p (r, r)G(-r, -r). (5) 

Here T pp and T s f are density and spin-fluctuation T-matrices, and T pp is the Cooper-pair 
fluctuation T-matrix. For U < and sufficiently low temperature, T pp provides the largest 
contribution to the self-energy; it is given by 

where Xpp( r 5 r ) = G(r, t)G(t, t). We have included all contributions to the FEA self-energy 
since this leads to a self-energy that is more accurate at high frequency (Deisz, et al. 1997). 

In Fig. 2, we show the superfluid density calculated for U = —4 and a density n = 0.75 
on a 64 x 64 and on a 128 x 128 lattice.^ The sharp upturn in D S (T) signals the transition 
to superconductivity; the finite size of the lattice is evident in the changes in D S (T) even for 
these large lattice sizes (Deisz et al. 1998a). As shown in the inset, the chemical potential 
becomes increasingly negative with increasing temperature as expected for a self-consistent 
conserving approximation that includes particle-particle 'ladder diagrams' (Serene 1989), as 
observed in quantum Monte Carlo calculations (Randeria et al. 1992), and in contrast to 
calculations in the 'ladder approximation' that are not self- consistent (Schmitt-Rink et al. 
1989). 

The single-particle spectral function A(k,e), provides a measure of the number of states 
accessible to an injected electron or hole. Roughly speaking, the product of A(k, e) and a 
Fermi function may be observed in angle-resolved photoemission experiments. The spectral 
function is related to the retarded propagator 

7r 7r (e — £k — Re £ (k, e)Y + (Im S K (k, e)) z 

We calculate A(k, e) from the retarded self-energy S R (k, e) which we obtain by analytic 
continuation from the imaginary frequency axis using Pade approximants (Vidberg and 
Serene, 1977). In contrast to quantum Monte Carlo data, our self-energies have no numerical 
noise of statistical origin, but do contain roundoff errors dependent on the criterion for 
convergence; we converged our self-energies so that the modulus of the largest change in E 
for any k point or Matsubara frequency from one iteration to the next is less than 1.0 x 10~ 7 
for 64 x 64 lattices and less than 1.0 x 10 -8 for 128 x 128 lattices. The density was held 
fixed at n = 0.75 to the same accuracy. 

In Fig. 3a, we compare single-particle spectral functions A(k, e) for several temperatures 
at one point inside but very near the Fermi surface k = (23,0)7r/32. The highest temper- 
ature shown (T = 0.20) for these calculations on a 128 x 128 lattice is well above any sign 
of superconductivity as shown in Fig. 2. At T = 0.20, A(ko,e) has a single rather broad 
peak.Q As the temperature is lowered in increments of 0.02, the peak in A(ko,e) shifts to 

4 We measure all energies and temperatures in units of the hopping matrix element t. 
5 The large width of the peak and the temperature dependence of the width is not consistent with a 
quasiparticle of a Landau Fermi liquid. We do not pursue this point here. 



5 




Figure 2: Supcrfluid density as a function of temperature for a 128 x 128 lattice (•) and 64 x 64 lattice (o). 
The arrows indicate the temperatures at which the excitation spectra are shown in panels 1 through 4 in Fig. 
1. The inset shows the self-consistent chemical potential as a function of temperature. The (temperature 
independent) Hartree contribution to the self-energy is included in the chemical potential. 

lower energies; at the lowest temperature shown (T = 0.12), spectral weight is suppressed 
at e = and it appears that the formation of a second peak for e > is incipient; cal- 
culations on 64 x 64 lattices at lower temperatures clearly show a second peak for e > 
(Hess et al 1996). Fig. 3b shows the real part of the denominator of the retarded Green's 
function, Re [G R (ko, e)] -1 = e — £k — ReS R (k ,e). The zero crossing of this expression, 
which corresponds to the energy of a quasiparticle excitation, occurs nearly at the peak 
in A(ko,e). With decreasing temperature the zero crossing shifts to lower energy. Struc- 
ture in Re [G R (ko, e)]" 1 at low energy also evolves with decreasing temperature. For the 
128 x 128 lattice, a perceptible superfluid density is only evident for the lowest temperature 
(see Fig. 2). To better understand the evolution of the peak for e > in A(k ,e) with 
decreasing temperature, we also consider Im T, R (k ,e). As seen in Fig. 3c, structure evolves 
in Im E fi (k , e) with decreasing temperature that is more complex than that which evolves 
in Re [G R (ko,e)]~ 1 . An examination of Fig. 3a shows that these changes do not lead to 
significant changes in the shape of the spectral function for the higher temperatures where 
the only apparent change in Re [G R (ko, e)]~ l is the shift of the zero crossing to lower energy. 
A more careful examination shows that while changes in the energy dependence of Im E are 
reflected in A(k , e), |Im E| is still rather large and changes of Re [G R (k Q , — Re£(k , e) 
dominate. Included in these are changes in the chemical potential required to hold the den- 
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sity fixed. Fig. 3d shows the spectral weight at e = as the temperature is decreased; a 
reduction of some 20% is evident as T c is approached and in the absence of any signature of 
sup er conduct i vity. 




Figure 3: Contributions to the spectral function for ko = (23,0)7r/32, a point near the Fermi surface 
shown for T = 0.20 (dash-dot), 0.18 (long-dash), 0.16 (dash), 0.14(dotted), and 0.12 (solid): (a) single- 
particle spectral weight, (b) Real part of the denominator of the fully rcnormalizcd retarded Green's function 
e — £k — ReE ii (k,£), and (c) Imaginary part of the self-energy. Also shown in (d) is the spectral weight 
at zero energy as a function of temperature. Comparison with the supcrfluid density and it's derivative 
in Fig. 2 suggests that the slight kink around T = 0.135 is another signature of the phase transition to 
superconductivity. 

The total DOS is a sum of A(k, e) over all momenta in the zone; our results suggest 
that the DOS will show a similar suppression of spectral weight. The suppression will 
differ quantitatively from that shown in Fig. 3d, which tracks the peak of what might be 
regarded as the 'coherent part' of the (non- Fermi liquid) quasiparticle excitation. Changes 
from incoherent parts of spectral functions for other k are included in the DOS. Depending 
in detail on how Im S(k, e = 0) changes with temperature for k near kp, these may lead to 
additional reduced contributions to the DOS as coherent quasiparticle excitations or possible 
gap structure in the low-energy excitation spectrum attempt to evolve. The evolution of 
these 'dimples' in the spectral function for k near kp as e — > would be analogous to those 
observed as a coherent state evolves in the Anderson lattice model (McQueen et al. 1993). 

These changes in single-particle excitation spectra are accompanied by the self-consistent 
formation of a sharp peak in T pp (q, uj m ) at q = and uj m = which evolves as an instability is 
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approached (see Fig. 3 of DHS). For larger lattices, T pp is more sharply peaked with a larger 
T pp (0, 0). Supposing that this peak leads to the dominant contribution to the self-energy,^ 
we find, 

T 

E(k,ie„) = — Y, G (- k + ^- £ n + ^m)T pp (q,uj m ) (8) 

m,q 

where i pp = T/N E| q |<£ T pp (q, 0), £ is an appropriate correlation length, and N is the number 
of sites. Inserting this result into Dyson's equation leads to a connection between G at (k, e n ) 
and the time reversal conjugate point (— k, — e n ) as expected for a pairing interaction. Taking 
k to be on the Fermi surface and using Dyson equations for G for two (k, e n ) points related 
by time reversal conjugation, we find for the single-particle spectral function, 

A(k F ,e) = -^y/d pp -£ 2 ; \e\ < 2^~. (9) 

This is a parabolic shaped spectral function centered on e = with a maximum value 
oc 1/ y^pp and with all spectral weight contained in the interval |e| < 2\Ji p ~ p . As t pp increases, 
the maximum is reduced and the size of the interval increases so that the sum of all spectral 
weight at any given k-point is unity. This is not the entire contribution to the self-energy, 
but this contribution is rapidly changing, leading to a suppression of spectral weight for any 
k-point at the Fermi surface. We note in passing that spectral weight is pushed away from 
the Fermi surface, and for k-points off of the Fermi surface, a second peak appears in the 
spectral function consistent with particle-hole coherent Bogoliubov-like excitations. While 
this model is too simplistic to correctly capture the essential features of the full calculation, 
it does illustrate how a superconducting transition in two dimensions that separates two 
disordered phases might occur in the FEA without giving rise to a finite gap but still leading 
to a suppression of low-energy spectral weight in the normal state. 

We return now to the excitation spectra shown in Fig. 1 for k along the (1, 0) direction 
in the zone, which shows that the quasiparticle excitations are still evolving as the super- 
conducting transition is approached, intermingled in a self-consistent way with the evolution 
of superconducting fluctuations. As the temperature is lowered, the suppression of spectral 
weight at zero energy is evident in the failure of a sharp peak to emerge in the narrow green 
region sandwiched between much sharper (red) quasiparticle and quasihole peaks slightly 
above and below zero energy. The largest value of the spectral functions occurs at the lowest 
temperature shown and is ~ 1.4. We have noted that the spectral weight at the point closest 
to the Fermi surface decreases strongly below T ~ 0.18 and falls by some 20% just above T c . 
So, relative to the sharp quasiparticle excitations above and below the Fermi surface, spectral 
weight at the peak of the quasiparticle excitation near the Fermi surface is suppressed by 
roughly 50%. At the lowest temperature, a finite superfluid density exists (see Fig. 2), finite 
but small spectral weight exists at zero energy, and hints of additional structure emerging 
from the region around zero-energy are evident. These "satellite" peaks suggest the coherent 
coupling of particle and hole excitations expected for a superconducting state, Bogoliubov 
excitations. 

6 A similar model was considered by Patton (1971). 
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While we have explicitly shown that spectral weight at low-energy is suppressed for this 
sizeable density, the quantitative result is less than observed in experiments for BSCCO 
which show a large reduction in spectral weight. The quantitative dependence of the spec- 
tral weight just above T c on the interaction strength and density is so far unknown. The 
dependence on the symmetry of the order parameter is also a subject of current research 
(Engelbrecht 1997). It is also largely unknown how Cooper-pair fluctuations interact with 
other possible mechanisms for the formation of the pseudogap that are not a consequence of 
superconductivity. Further work seeking to elucidate these issues is underway. 

We have presented fully self-consistent calculations of the temperature dependent su- 
perfluid density together with the single-particle excitation spectrum for a strong coupling 
superconductor with an attractive interaction of half the bare bandwidth in two dimensions 
on large 128 x 128 lattices. We have taken the electronic density to be roughly that expected 
from electronic structure calculations. Our calculations show the evolution of quasiparticle 
excitations self-consistently intermingled with superconducting fluctuations, and explicitly 
demonstrate a significant suppression of spectral weight at low energy as the superconducting 
transition is approached from above. For temperatures just below the transition temperature, 
a gap is not evident. A simplistic self-consistent model was presented to illustrate how the 
FEA might produce a phase transition between disordered states in two dimensions without 
producing a finite gap in the single particle excitation spectrum. 
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